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Abstract 

The instability of a Fermi-liquid drop with respect to bulk density distortions 
is considered. It is shown that the presence of the surface strongly reduces 
the growth rate of the bulk instability of the finite Fermi-liquid drop because 
of the anomalous dispersion term in the dispersion relation. The instability 
growth rate is reduced due to the Fermi surface distortions and the relaxation 
processes. The dependence of the bulk instability on the multipolarity of the 
particle density fluctuations is demonstrated for two nuclei Ca and 208 Pb. 
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I. INTRODUCTION 



In the vicinity of the equilibrium state of the nuclear Fermi-liquid drop the stiffness co- 
efficients are positive and the system is stable with respect to particle density and surface 
distortions. With decreasing bulk density or increasing internal excitation energy (temper- 
ature) the liquid drop reaches the regions of mechanical or thermodynamical instabilities 
with respect to small particle density and shape fluctuations and to separation into liquid 
and gas phases. The process of the development of instability is a complicated one. We will 
discuss some aspects of the process. In particular, we will study the influence of a sharp 
liquid-drop boundary on the instability with respect to small particle density fluctuations. 
In actual nuclear processes (heavy-ion reactions, nuclear fission etc.) nuclear matter is not 
static, and consequently the development of instability depends not only on the equation of 
state, but also on the dynamical effects such as the dynamical Fermi-surface distortion or 
the relaxation processes. We will take into account these aspects in studying the stability 
of the Fermi-liquid drop in both regimes of the first- and zero sound modes. 



II. BULK INSTABILITY OF THE FERMI-LIQUID DROP 

Let us consider small density fluctuations 5p(r, t) starting from the nuclear fluid dynamic 
approach [PlEf. The linearized equation of motion reads (see Ref. [[|), 



Q2 foj^ 

m-^Sp = Vp eq V— + V„V„P^, (1) 



where p eq is the equilibrium density, E is the total energy and the pressure tensor, P' , 
represents the deviation of the pressure from its isotropic part due to the Fermi surface 
distortions. 

The variational derivative 5E/Sp in Eq. (P implies a linearization with respect to the 
density variation Sp: 



We point out that the first term on the r.h.s. of Eq. (0) does not enter Eq. ([!]) because of 



the the equilibrium condition (5E/5p) eq = Xp = const, where is the chemical potential. 
The operator L can be derived from the equation of state E = E[p]. We will use the 
extended Thomas-Fermi approximation for the internal kinetic energy |J and the Skyrme- 
type forces for the interparticle interaction [|J. In the special case of a spin saturated and 
charge conjugated nucleus, and neglecting spin-orbit and Coulomb effects, the equation of 
state reads 
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The effective forces used in Eq. ([!]) leads to an overestimate of the incompressibility co- 
efficient. This is a well-known feature of Skyrme forces which can be overcome by taking 
non- integer powers of p in the potential energy density in Eq. ([3]). For our purposes we shall, 
however, be content with the form (|3|). To make quantitative estimates of the finite size ef- 
fects on the bulk instability of the liquid drop, we will assume a sharp surface behaviour of 
p e q(r) having a bulk density po and an equilibrium radius Rq. Taking into account Eqs. @ 
and (H), the operator L[p eq ] is then reduced to the following form 
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and K is the incompressibility coefficient 
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The Landau parameters Fi are given by 
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where 



m m on 
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and is the Fermi energy. 

The pressure tensor P' can be expressed through the displacement field x(r, t) |§ . 
Assuming also Sp ~ e~ ia; *, the pressure tensor P' is given by [||] 



where r is the relaxation time and we used the symbol 

2 . , 

= V^Xm + V M x„ - 2^m V aXa (8) 

for this combination of gradients of the Fourier transform \ v of the displacement field x(r, t). 
The equilibrium pressure of a Fermi gas, P eq , in Eq. (0), is given by 

1 f dp 



P - = ^/?2^ p2/ - (r ' P) " P0 ^ /5m ' 



3m J (2nhy 

where f eq (r, p) is the equilibrium phase-space distribution function and pp is the Fermi 
momentum. We point out that Eq. (0) is valid for arbitrary relaxation time r and thus 
describes both the zero- and the first-sound limit as well as the intermediate case. 

Taking into account the continuity equation and Eqs. (f|), (^) and (p]), the equation of 
motion (|1|) can be reduced in the nuclear interior to the following form (we consider the 
isoscalar mode): 

-mu 2 5p= { X -K-\ -^r (P eq /po)) V 2 5p -2(J3 + t a p ) V 2 V 2 5p. (9) 



,9 3 1 — iujt 
The solution of Eq. (|9|) for a fixed multipolarity L is given by 

Sp(r,t) = p j L (qr) Y LM (6,<f)) a LM (t), (10) 

where q is the wave number and ola/(£) is the amplitude of the density oscillations. We 
will distinguish between stable and unstable regimes of density fluctuations. In the case of 
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a stable mode at K > 0, a solution of Eq. (|9|) of the form (|T0|) has the following dispersion 
relation 

2 2 2- 70*0 2,4 

u) =u q —iuj q +n s q . (11) 

m 

Here, u is the sound velocity 

u 2 = u\ + «„, (12) 
where Mi is the velocity of the first sound 

uj = —K, (13) 

j(uj) is the viscosity coefficient 

= (14) 
3 v 1 — iujt j p eq 

and 

= ^Im(— — — )-^ SL , k s = — (/? + t s p ). (15) 
3 v l — iuot / m p eq m 

The quantities and 7(0;) appear due to the Fermi-surface distortion effect. The dispersion 
relation (|ll]) determines both the real and the imaginary part of the eigenfrequency uj. 

The equation of motion @ has to be augmented by the boundary condition. This is 
given by a condition of the balance of the surface pressure 5P sur f with the volume sound 
pressure SP soun d on a free surface of the liquid drop, see Refs. It reads 

Wpo jL(qRo) = ^1 (L - 1} (L + 2) ' ^a^U,- (16) 

Let us consider now the volume instability regime, K < 0, and introduce a growth rate 
r = — ioj (r is real, T > 0), see Ref. 0. Using Eq. (|TT|), one obtains 

T 2 = \ Ul \ 2 q 2 -ar)q 2 -K s q\ (17) 

where 



Equation (|17D is valid for arbitrary relaxation time r. From it one can obtain the leading 
order terms in the different limits mentioned above. 

(i) Frequent collision regime: r — > 0. 
The contribution from the dynamic distortion of the Fermi surface, k v , can be neglected in 



this case and we have from Eqs. (|ITD and (14]), 



r 2 = \ Ul \ 2 q 2 -T (7/m)g 2 - K s q\ (19) 

where 7 = (8/15) cft is the viscosity coefficient. In the case of small viscosity coefficient 7, 
one has from Eq. fllTf ) 
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T 2 ^ \ Ul \ 2 q 2 - k s q 4 - — q 2 J\ Ul \ 2 q 2 - K s q*. (20) 
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The amplitude of the density oscillations, 5pl(t, t), grows exponentially if T > 0. Ex- 
pression ( |20"D determines two characteristic values of the wave number q, namely, q m ax where 
the growth rate reaches a maximum of T max , and q crit where T goes to zero, i.e., (see also 

El), 

r = T max at q = q max < q crit , and T = at q = q crit . (21) 



\ u |2 

= and q 2 crit = L -^-, at u\ < 0. (22) 

Q — Qmax 



The values of q max and q cr u are obtained from, see Eq. (|T9"|) 

ar 

Thus, the critical wave number q crit does not depend on the viscosity. However, the presence 
of viscosity reduces the instability, see also Fig. 1 below. 

(ii) Rare collision regime: t — > 00. 
In this case, we have from Eqs. (|I7D , (|15D and (0) 

T 2 = \u 1 \ 2 q 2 -K' v q 2 -K s q\ (23) 

where 
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The critical value q crit and the value q max are given by 
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Thus, the distortion of the Fermi-surface leads to a decrease of the critical value qcrit, i-e., 
the Fermi-liquid drop becomes more stable with respect to the volume density fluctuations 
due to the dynamic Fermi-surface distortion effects. 

III. NUMERICAL RESULTS AND DISCUSSION 

In Fig. 1 we have plotted the instability growth rate V as obtained from Eq. flTTl). The 
calculation was performed for the Skyrme force SIII. The relaxation time was taken in the 
form r = ha/T 2 with a = 9.2 Me V and a = 2.6 Me V [0 and the bulk density po wa s 



taken as p = 0.3 p sat , where p sat is the saturated density p sat = 0.1453 fm~ 3 . We show also 
the result for the nonviscous infinite nuclear matter and the nonviscous finite liquid drop 
neglecting Fermi surface distortion effects. In a finite system, the non-monotony behaviour 
of the instability growth rate as a function of the wave number q is due to the anomalous 
dispersion term in Eq. ( |TT1) created by the gradient terms in the equation of state. We 
point out that the finite system becomes more stable with respect to short-wave-length 
density fluctuations at q > q max . We can also see that the presence of viscosity decreases the 
instability. The strong decrease of instability in a Fermi liquid drop (FLD), when compared 
with the corresponding result for the usual liquid drop (LD), is because of the Fermi surface 
distortion effects. In Fig. 2, this peculiarity of the FLD can be seen in a transparent way 
for both the infinite nuclear matter and the finite Fermi liquid drop. 

For a saturated nuclear liquid one has for the force parameters t < 0, t 3 > and t s > 0. 
Thus, the critical value q cr u, Eq. (p2|), increases with decreasing bulk density po at u\ < 0, 
see also Eq. (|5]). The existence of the critical wave number q crit for an unstable mode is a 
feature of the finite system. The growth rate T depends on the multipolarity L of the nuclear 



density distortion and on the position of the eigenvalue, qL, in the interval of q = -j- q cr it 
||10||. For a given i? 0) the value of qi increases with L for L > 2 because of the boundary 



condition (jig), see Table 1. That means that if < g max the instability increases with L 
and the nucleus becomes more unstable with respect to an internal clusterization to small 
pieces (high multipolarity regime) rather than to binary fission (low multipolarity regime). 
In contrast, the binary fission is preferable if q max < qi < q C rit- 

We give in Table 1 the values qhj^F for two nuclei, 208 Pb and 40 Ca as obtained from 
Eq. (|i~6|) . The calculations were performed with the surface tension parameter Airr^a = 
17.2 MeV. We point out that the value of q max is given here by q ma x/kF = 0.69. In Fig. 
3 we have plotted the instability growth rate at T = 6 MeV and a = 9.2 MeV as function 
of the multipolarity L of the particle density fluctuations for two nuclei 208 Pb and 40 C a. 
As is seen from Fig. 3, the lowest values of L < 3 give the contribution to the instability 
growth rate Y for the nucleus 40 C a. Thus, the nucleus 40 C a is unstable with respect to the 
fission under the conditions considered above. In contrast, the instability growth rate of the 
nucleus 208 Pb includes the higher multipolarity L < 8 and this nucleus has to be unstable 
with respect to multifragmentation. 



IV. SUMMARY AND CONCLUSION 

Starting from the fluid dynamic equation of motion for the Fermi liquid drop with a 
sharp surface, we have derived the dispersion relations (|TT|) and (|P7|) for both the stable and 
the unstable regime. The dispersion relations are influenced strongly by the Fermi-surface 
distortion effect and the anomalous dispersion caused by the finiteness of the system. The 
presence of the Fermi surface distortion enhances the stiffness coefficient for a stable mode 
and reduces the instability growth rate for an unstable one. 

We have shown that the instability growth rate in an unstable finite system is a non- 
monotony function of the wave number q because of the anomalous dispersion term. This is 
in contrast with the infinite nuclear matter case where the instability growth rate increases 
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with q. The non-monotony behaviour of the instability growth rate T(q) in a finite Fermi 
liquid drop is accompanied with two characteristic wave numbers q max and q cr it, see Eqs. 
(pl|) and (p2|). The distortion of the Fermi-surface leads to a decrease of the critical value 
q cr it- The decay mode of an unstable Fermi liquid drop depends on the location of the eigen 
wave number q on the slope of the curve T(q). The Fermi liquid drop is more unstable with 
respect to multifragmentation if q < q max and the binary fission is preferable if q > q ma x- 
This is because the eigen wave number qi, derived from the secular equation (|16|) , increases 
with the multipolarity L of the particle density fluctuations. As an example, we have 
demonstrated this phenomenon in the case of hot nuclei m Ca and 208 Pb. The nucleus m Ca 
is more unstable with respect to the short wave fluctuations and prefers to decay into the 
binary fission channel. The multifragmentation channel is preferable at the development of 
the instability in the heavy nucleus 208 Pb. 
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TABLES 

TABLE I. Values of qi/kp, as obtained from Eq. (16) for the surface tension parameter a 
derived from \ -Krla = 17.2 MeV and T = lMeV, for two nuclei 208 Pb and 40 Ca. 



L 


2 


3 


4 


5 


6 


7 


8 


40 Ca 


0.986 


1.225 


1.457 


1.684 


1.907 


2.128 


2.347 


208p b 


0.569 


0.707 


0.841 


0.972 


1.101 


1.228 


1.355 
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FIGURES 

FIG. 1. The dependence of the instability growth rate T on the wave number q. The calculations 
were performed for Skyrme force SIII, temperature T = 6 MeV and density po = x p sa t with 
x = 0.3. The solid curves are for the Fermi liquid drop from Eq. (17); the values a = 9.2 MeV 
and a = 2.6 MeV of the relaxation time parameters are shown as labels to the curves. The dashed 
lines are the results for the nonviscous liquid without the Fermi surface distortion effects: curve 
(1) is the result for an infinite matter and curve (2) is for a finite liquid drop. 

FIG. 2. The dependence of the instability growth rate T on the wave number q for the infinite 
matter (dashed lines) and the finite system. The curves FLD and (2) are for the Fermi liquid; the 
curves LD and (1) are for the usual liquid, i.e., neglecting the Fermi surface distortion effects. All 
calculations were performed with a = and T = 1 MeV and with the force parameters and po as 
in Fig. 1. 

FIG. 3. The dependence of the instability growth rate T on the multipolarity L of the particle 
density fluctuations for two nuclei 208 Pb and m Ca. The calculations were performed using the 
FLD results of Fig. 1 at a = 9.2 MeV and T = 1 MeV and the surface tension parameter a derived 
as in Table 1. 
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